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Abstract The possibility to use the ultrashort ion bunches produced by 
circularly polarized laser pulses to drive a source of fusion neutrons with 
sub-optical cycle duration is discussed. A two-side irradiation of a deuter- 
ated thin foil target produces two countermoving ion bunches, whose col- 
lision produces an ultrashort neutron burst. Using particle-in-cell simula- 
tions and analytical modeling, it is evaluated that, for intensities of a few 
10^^ W cm^^, more than 10'^ neutrons per Joule may be produced within 
a time shorter than one femtosecond. Another scheme based on a layered 
deuterium-tritium target is outlined. 
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1 Introduction 

With the advent of short-pulse laser systems yielding multi-terawatt power, 
laser-driven nuclear physics has emerged as a very active area of research 
with applications such as radioactive isotope production and nuclear trans- 
mutation of elements. In particular, the emission of neutrons from fusion 
reactions has been observed in several experiments with different pulse pa- 
rameters and target types, including solid targets heavy water droplets 
[3], deuterium clusters and underdense plasmas or gas jets 6 . In these 
experiments, the number of neutrons produced per Joule of the laser pulse 
energy is usually in the 10* — 10^ J^^ range for terawatt, femtosecond lasers, 
and higher for large petawatt, picosecond lasers (see table I of Ref.j3 for 
a partial summary). While the large neutron fluxes produced by petawatt 
pulses are of importance for material damage studies relevant to thermonu- 
clear fusion research [7], neutron bursts produced by "table-top terawatt" 
(T'^) lasers at high repetition rate may provide compact, pulsed neutron 
sources for radiography and other applications. 

The interpretation of the above experiments suggests that, as a general 
issue, the route to neutron production starts from the heating of target elec- 
trons up to high energies; in turn, the electron currents produce space-charge 
fields driving ion acceleration up to MeV energies (via various mechanisms 
such as sheath acceleration in solid targets, or Coulomb explosion in clus- 
ters and underdense plasma channels); finally, collisions between ions lead 
to fusion reactions and neutron emission. The duration of the latter could 
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not be measured so far in experiments; however, based on the physical pic- 
ture which is inferred from the experiments, we expect the neutron pulse 
duration not to be shorter than the laser pulse duration, i.e., to be typically 
in the 0.1 1.0 ps range. 

In this paper, we study a new approach to neutron production by ul- 
trashort laser pulses, which could provide a source of fusion neutrons with 
sub-optical cycle duration. This represents another possible route to the 
ultrafast control and imaging of nuclear reactions by superintensc fields |S] 
and may be used to study phenomena such as nuclear spin-mixing oscilla- 
tions whose period was estimated to be ^ 1 fs 0. 

Our approach is based on the prediction that the irradiation of a solid 
targets with circularly polarized, intense short pulses leads to the prompt 
acceleration of high-density, short duration ion bunches jl()| . Such bunches 
may be used efficiently as projectiles for beam fusion, since for deuterium the 
energy in the center-of-mass system may approach the Gamow value {£g « 
1 MeVmr/mp, where rrir is the reduced mass) for which the cross sections 
of deuterium-deuterium (D-D) or deuterium-tritium (D-T) reactions have a 
maximum. With an appropriate target scheme, fusion reactions and related 
neutron emission may last just for a time of the order of the ion bunch 
duration, which may be less than one optical cycle. 
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2 The colliding bunches scheme 

To find a scheme based on the D + D ^ ■^He + n (2.45 MeV) reaction aiming 
at the shortest achievable duration, we consider a symmetrical, double-sided 
irradiation of a thin foil target. In this scheme, two colliding ion bunches are 
generated. Thus, if the two bunches are properly timed, for a given value 
of the laser intensity the energy in the center-of-mass is maximized while 
the duration of neutron emission is minimized. This experimental geometry 
is similar to the one of Ref.|llj. where a "laser-confined" thermonuclear 
fusion approach was proposed. In that scheme, a thin deuterium-tritium 
foil is compressed, heated and confined by double-side irradiation using 
relatively long pulses. In our scheme, to produce two single bunches the 
laser pulse duration must be much shorter (a few cycles) than in R.ef.|ll|. 
Our approach is based on a non-thermal beam fusion concept with emphasis 
on the ultrashort duration of the neutron emission, and without concerns 
of target stability. 

We now discuss the target and pulse requirements based on the features 
of ion bunch generation studied in Ref . |l(Jj . The bunch ion velocity spectrum 
extends up to a maximum velocity given by 

!^ ^ 2 /— — — fli ~ 047 /^^fli (1) 

C V ^ "^P \ A Tie ' 

Tie is the initial electron density, — meCt;£/47re^ is the cut-off density for 
the laser frequency ujl, and ql — O.SSy^/A^/lO^^ W cm^^/^m^ is the dimen- 
sionless laser amplitude, being / and A ~ i'Kc/ujL the intensity and the wave- 
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length of the laser pulse, respectively. At a given value of a^, the lower the 
higher Collisions in the center-of-mass system at the Gamow energy for 
the D-D reaction will begin when 2vm = vq = \J 2£g / mp ~ 0.0458c. The 
ion bunch acceleration occurs over a typical time ~ ITj^i^AjZYl'^aL'^ 
where is the duration of a laser cycle. Thus, ultraintense few-cycle pulses, 
which are now at the frontier of current research |12| . are best suited for ion 
bunch acceleration. 

The number of accelerated ions per unit surface is ~ riioh, where n^o 
is the background ion density and Ig is the evanescence length of the pon- 
deromotive force inside the plasma. At the end of the acceleration stage 
(i = Ti), a sharp peak of the ion density is formed at a distance ~ Ig from 
the original target surface; the peak density is Ub = kuq^ where k ~ 10 is 
found from simulations, and the bunch width is thus ~ Is/k. One ex- 
pects Is — c/ojp = (A/27r)^ Tic/ne, where ujp — ^ Annee^ / nie is the plasma 
frequency. Thus, /f, <C A holds when rig ^ ?ic, and the bunch duration 
Tf, ~ h/vm and can be less than = A/c. 

In the colliding bunches scheme the duration of the neutron burst, will 
be of the order of Tb/2, i.e. in the sub-femtosecond regime. The optimal 
target thickness is close to the value £ ~ 2ls, in order to let the two bunches 
collide at the end of the acceleration stage. Since a very thin foil target 
is needed, deuterated plastic, e.g. (CD2)n, looks more suitable than pure 
(cryogenic) deuterium as a target material; however, the much higher value 
of rie requires higher pulse intensities. In the following, both pure solid D 
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(rii = rie ~ no = 6 X 10^^ cm^'^ i.e. Ue = AOric for A = 0.8 /im) and "plastic" 
targets (with same = no but Ue = 250nc) are discussed. 



3 Simulation results 

We now analyze the double-sided irradiation of a thin foil using a particle- 
in-cell (PIC) simulation. As shown in Ref. Ji)_ a one-dimensional (ID) model 
is adequate to keep the essential features of ion bunch generation. This al- 
lows us to use a very fine spatial and temporal resolution, which is actually 
necessary to resolve the bunch formation properly, since the latter is char- 
acterized by sharp density gradients and very short formation times. In our 
simulations we use 2000 grid cells per wavelength and 2000 particles per 
cell at the initial time. In addition, to compute the neutron yield from the 
simulation data (as shown below) properly, the ion distribution function has 
to be reconstructed over a phase space grid once every a few timesteps. 

Figurenshows results for a "D" foil simulation. The front and rear laser 
pulses both have peak amplitude — 2.5 and duration tl — 5Tl (FWHM), 
corresponding to an intensity of 1.3 x 10^^ W cm^^, and a duration of 13 fs if 
A = 0.8 /im. The foil target has initial density rie — n^o = 40nc and thickness 
£ — 0.05A. The ion density reaches a peak value close to 1390nc ~ SSn^o, 
i.e. K ~ 17, with rise and fall times of sa O.ITl- The maximum ion velocity 
Vm — 10^'^c. The results of a "CD2" simulation {rie/ric — 250, gl = 8, 
£ = 0.02A and all the other parameters equal to the case of Fig^ are 
qualitatively similar: k ~ 10 and Vm — 0.01c are found. 
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We now compute the neutron burst intensity and duration from the 
simulation data. The reactivity for a fusion reaction is given by 



where <y{v) and g{v) are the cross section and the distribution function 
(normahzed to unity), respectively, in terms of the relative velocity v — 
|vi — V2I, being vi.2 the velocities of the two ion species. The number of 
fusion reactions per unit volume and time is given by i? = nin2 {o'v)/{l+Si2) 
|13| . being ni_2 the number densities of the two ion species. For ions of 
the same species, as in the D-D reaction, ni — n2 — n. A convenient 
parameterization of the cross section is given by |13j 



where £ = is the center-of-mass energy, and S is the astrophysical 

factor that is a slowly varying function of £ and hence will be taken as a 
constant 5*0 = 5.4 x lO^^s i^^y j^ j^ 

From the PIC simulation data we obtain g{v) at any grid point by com- 
puting the ion velocity distribution f{vi) on a phase space grid, from which 
g{v) is obtained by convolution. We thus compute the reactivity and, from 
the knowledge of the ion density, the number of fusion reactions per unit 
density and time at each grid point. Finally, we obtain the total rate of 
fusion events and the overall number of neutrons produced by integrating 
R over space. 




(2) 



a = 



c ^ ' 



(3) 
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The results are shown in FigI21for both the "D" and "CD2" target cases. 
The same number density of deuterium ions has been assumed. In both 
cases 5 a neutron burst is generated with a duration of about 0.7 fs (FWHM) 
and a yield of ^ lO^ncutrons cm^^. Since the pulse duration is ~ 15 fs, the 
number of neutrons produced per Joule of the pulse energy is ^ 10^ J^^ for 
the D case with QL = 2.5, - 10^ CD2 case (fli = 8); these numbers are 
roughly between one and two orders of magnitude lower than that inferred 
from experiments with T'^ systems; however, for the present scheme the 
expected duration is likely to be much shorter, and the emission rate may 
be comparable or even higher. 

The relatively long-lasting tail in the fusion rate is due to the fact that 
the two bunches, after crossing each other, propagate in a low-density shelf 
of ions which originate from the layer of charge depletion at the surface and 
have been accelerated to lower energies than the bunch ions (see Fig^and 
E,ef.|l()|). At higher intensities the ions in the shelf have enough energy to 
sustain a significant rate of fusion reactions during the propagation of the 
bunch, leading to a higher total yield but also to a longer duration (a few 
femtoseconds) of the neutron emission. 

4 Analytical scalings 

To check our numerical results, and to provide an approximate scaling of the 
neutron yield versus the laser intensity, we analytically evaluate the fusion 
rate in the colliding-bunches scheme. From Ref.pO] we infer that it is worth 
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to consider for the fon bunch vefocity distributfon f{vi) either a "flat-top" 
velocity distribution fp = l/w„i for < Vi < Vm and zero elsewhere, and 
a monochromatic distribution ~ 5{vi — Vm). In both cases the relative 
velocity distribution g{v) is easily obtained analytically by convolution. The 
averaged reactivity may thus be written as 

450 f 6-"^/" , , , _ 45'o 



iav) 



g{v)dv= ^A^(C), (4) 



mpVG J v/vg rripVG 
where C = wc/wm, 

f 2C [Ei(C/2) - Ei(C) - E2(C/2) + E2(C)] Af = If) , _ 
M{0 = \ (5) 

[(C/2)e-^/^ if = fD), 

and En{x) — {e^^'^ /t^)dt is the exponential integral function of order n. 

The rate of neutrons produced per unit volume is given by i? = (n| /2){av), 
where n,; — kuq. To obtain the number N of neutrons produced per unit sur- 
face we multiply R by the spatial extension of the neutron burst, ~ Is/n, 
and by the burst duration ~ rb/2 ~ Is/lnvm- Hence, the total neutron yield 
does not depend on the compression ratio k, which only affects the bunch 
width and the burst duration. We obtain 

N li-iLl^ J^A^(c) ^ N,CM{0 . (6) 

2 2VG Vm rUpVG 

The dependence of N on the laser intensity is shown in Fig|31for both choices 
of f{v). Posing Is ~ c/ojp, we obtain iVo ~ 2 x 10^ for "D", A^o - 3 x lO"^ for 
"CD2" targets. These values underestimate the neutron yield observed in 
simulations; this is likely to be due to the fact that the effective screening 
length Is is actually larger than c/up and/or the active volume for fusion 
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reactions is wider than the bunch, since, as observed above, fusion reactions 
occur also in the low-density shelf. This explanation is also supported by 
noting that the burst duration is also underestimated when taking Is ~ c/ujp: 
for "D" targets, c/{2ujpKv„i) ~ 0.2 fs while a duration of 0.8 fs is observed 
in Figia 

5 Single bunch scheme 

Using ultrashort ion bunches, it may be also possible to obtain a fetmosec- 
ond source of neutrons or other fusion products from hcteronuclear reactions 
using a single short pulse impinging on a layered target. As an example we 
consider neutron production in a target with a thin surface layer of deu- 
terium for ion acceleration and an inner tritium layer as an immobile target, 
using the D + T ^ a + n (14 MeV) reaction, hi the relevant energy range^ 
the astrophysical factor for the D-T reaction is 5*0 = 1.2 x 10^^" keV cm^, 
almost 200 times larger than for D-D. This effect compensates the loss in 
the center-of-mass energy with respect to the colliding bunches scheme at 
the same laser intensity. If the deuterium layer has a thickness ~ Is, and 
the tritium layer is not thicker than the ion bunch, the neutron burst dura- 
tion might be limited to ~ It/vm- An analytical estimate, analogous to that 
^ The D-T cross section has a broad maximum around 64 keV; however, this 
low-energy range (where S{£) varies strongly with energy) is not considered here, 
since the laser-plasma interaction regime becomes very coUisional and must be 
further investigated. 
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leading to Eq.ljHl, leads to the following result 



n [To 



( 



Sq 



) 



DT 



■AiO ^ 1.3 X 10" cm-^K-\AiO , (7) 



where -4(C) = aCEi(aC) or a(e for a flat-top or monochromatic ion 



assumed. If the tritium layer It is thicker than 4, both the yield and the 
duration increase by a factor ~ Ix/h- The function QA{C) has a maximum 
in the range of intensities 10^° ~ 10^^ W cm~^. 

6 Discussion and conclusions 

We briefly discuss the experimental feasibility of the proposed scheme. The 
required laser pulse duration and intensity are within present-day, or at 
least near-term capabilities of table-top systems. Quenching of prepulse, 
pulse synchronization and need of circular polarization as well as thin foil 
target manifacturing are likely to be demanding tasks but not out of reach. 
The measurement of the duration of a neutron burst with sub-femtosecond 
resolution is a challenging issue. We argue that it may be possible to perform 
an indirect measurement based on the interaction of the neutrons with a 
secondary target, where charged particles are produced; this process may 
be resolved in time by using attosecond spectroscopy techniques jl4j. 

In conclusion, we investigated the production of fusion neutrons using 
ion bunches accelerated in solid deuterated targets by circularly polarized 
laser pulses. Our results, in connection with the continuous progress in pro- 
ducing ultrashort, superintense laser pulses at a high repetition rate, may 



distribution, respectively, a 



■\/l25/24, and no ~ ut — has been 
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open a perspective for a neutron source with sub-optical cycle duration and 
sufficient brightness for applications. 
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Fig. 1 (Color online.) The ion density rii (top row, solid line) and the phase space 
distribution f(x,Vx) (bottom row, arbitrary units) at different times (labels) from 
an ID PIC simulation of two-side irradiation of a thin foil. The initial density 
profile is also shown (dashed line). Run parameters are a_L = 2.5, no/iic = 40, and 
I = 0.05A. 
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Fig. 2 Numerical evaluation of the rate (solid line) and the total number (dotted 
line) of neutrons produced per unit surface for the "D" simulation of FigQ (left 
panel) and for the "CD2" simulation for which = 250nc, a_L = 8. 
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Fig. 3 Analytical scaling of the number of neutrons (per unit surface) produced 
during the ultrashort burst as a function of laser intensity for a A = 0.8 /xm 
pulse, for a "D" target with rie = 40nc and a "CD2" target with rie = 250nc. 
Solid and dotted curves correspond to the choice of "flat-top" and monochromatic 
distribution functions, respectively. 



